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Abstract 

This is the final fifth paper in our series of five in which we test the Master Constraint 
Programme for solving the Hamiltonian constraint in Loop Quantum Gravity. Here we 
consider interacting quantum field theories, specifically we consider the non - Abelean Gauss 
constraints of Einstein - Yang - Mills theory and 2+1 gravity. Interestingly, while Yang - 
Mills theory in 4D is not yet rigorously defined as an ordinary (Wightman) quantum field 
theory on Minkowski space, in background independent quantum field theories such as Loop 
Quantum Gravity (LQG) this might become possible by working in a new, background 
independent representation. 

While for the Gauss constraint the Master constraint can be solved explicitly, for the 
2+1 theory we are only able to rigorously define the Master Constraint Operator. We 
show that the, by other methods known, physical Hilbert is contained in the kernel of the 
Master Constraint, however, to systematically derive it by only using spectral methods is as 
complicated as for 3+1 gravity and we therefore leave the complete analysis for 3+1 gravity. 
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1 Introduction 

We continuue the test of the Master Constraint Programme for Loop Quantum Gravity 
(LQG) 13 El which we started in the companion papers 0E1IH IS]. The Master Constraint 
Programme is a new idea to improve on the current situation with the Hamiltonian constraint 
operator for LQG 9 . In short, progress on the solution of the Hamiltonian constraint has been 
slow because of a technical reason: the Hamiltonian constraints themselves are not spatially 
diffeomorphism invariant. This means that one cannot first solve the spatial diffeomorphism 
constraints and then the Hamiltonian constraints because the latter do not preserve the space 
of solutions to the spatial diffeomorphism constraint ^U]- On the other hand, the space of 
solutions to the spatial diffeomorphism constraint ^0] is relatively easy to construct starting 
from the spatially diffeomorphism invariant representations on which LQG is based jllj which 
are therefore very natural to use and, moreover, essentially unique. Therefore one would really 
like to keep these structures. The Master Constraint Programme removes that technical obstacle 
by replacing the Hamiltonian constraints by a single Master Constraint which is a spatially 
diffeomorphism invariant integral of squares of the individual Hamiltonian constraints which 
encodes all the necessary information about the constraint surface and the associated invariants. 
See e.g. 0] for a full discussion of these issues. Notice that the idea of squaring constraints is 
not new, see e.g. |13j . however, our concrete implementation is new and also the Direct Integral 
Decomposition (DID) method for solving them, see [2] for all the details. 

The Master Constraint for four dimensional General Relativity will appear in |15j but before 
we test its semiclassical limit, e.g. using the methods of ^§] E] and try to solve it by DID 
methods we want to test the programme in the series of papers [UGDIUIH]. m the previous papers 
we focussed on finite dimensional systems of various degrees of complexity and free field theories. 
In this paper we consider interacting quantum field theories, in particular in four dimensions and 
2+1 gravity which is interacting before one solves the constraint (the reduced phase space is of 
course finite dimensional). This might be a surprise given the fact that no interacting Wightman 
quantum field theories have been constructed so far in four dimensions. The resolution of the 
puzzle is that instead of working on a Minkowski background we couple gravity to matter and 
therefore arrive at a background independent quantum field theory. That these display a better 
ultraviolet behaviour than background dependent quantum field theories, at least for a large 
class of operators was already shown in [5] and we will see that the same mechanism is at work 
here and makes the associated Master Constraint Operators well defined. This does not mean, 
however, that we are done and have a consistent quantum theory of both matter and geometry 
because only 1. after having solved the Master Constraint for all constraints (here we only 
consider the matter Gauss constraints but not, in particular, the Hamiltonian constraints) of 
gravity coupled to the standard model and 2. after having shown that the solution theory is 
consistent with experiment (reduces to QFT on (curved) backgrounds in the semiclassical limit 
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of low geometry fluctuations) do we have a viable candidate theory. It is precisely the purpose of 
the Master Constraint Programme to complete this task and the present series of papers serves 
the purpose to make sure that the Master Constraint Programme reproduces the established 
results in solvable cases. 

2 Infinite Number of Non-Abelean First Class Constraints Non 
— Polynomial in the Momenta with Structure Constants 

In a previous article [H] we discussed a background dependent quantum field theory, specifi- 
cally free Maxwell theory. Even in this free theory the definition of the Master Constraint had 
to involve a non-trivial integral kernel C which depended on the background metric in oder 
that the Master Constraint Operator was densely defined. In this section we discuss two back- 
ground independent theories, namely Einstein - Yang - Mills theory and pure Einstein gravity. 
Notice that these theories are no longer free, both the Yang - Mills and gravitational degrees 
of freedom are self-interacting and interacting with each other. One therefore expects severe 
ultraviolet divergence problems, the more as we are not allowed to use a background metric in 
order to regulate the Master Constraint Operator. Nevertheless it turns out that precisely 
because of the background independence the ultraviolet problems can be overcome, 
background independent theories regulate themselves! The way this works is identical 
to the mechanism discovered in jHj for the definition of various Hamiltonian constraints in Loop 
Quantum Gravity (LQG) (see [7J|S] for an introduction). Hence we will construct a correspond- 
ing Master Constraint Operator for the Gauss constraints of Einstein - Yang - Mills theory 
using LQG techniques. 

In what follows we consider canonical Yang - Mills fields for a compact gauge group G 
with canonicaly conjugate pair (A£,E_j). Here a,b,.. = 1,2,3 are spatial tensor indices while 
J, K, .. = 1, ..,dim(G) are Lie(G) indices. Likewise for the gravitational sector we consider the 
canonical pair (Aa,Ej) where j,k,.. = 1,2,3 are su(2) indices. From the kinematical point of 
view the gravitational phase space is identical with that for a Yang - Mills theory with gauge 
group SU(2). We assume both principal bundles to be trivial for simplicity. The non- vanishing 
Poisson brackets are 

m(*)>4*(v)} = 9 2 S%6f6(x,y) 

{E*(x),A k b (y)} = K5 a b 5*5(x,y) (2.1) 

where g 2 and k = 8irG Newton denote the Yang - Mills and gravitational coupling constant 
respectively. We are using units in which both connections have dimension cm -1 while the Yang 
- Mills and gravitational electric fields respectively have dimension cm' 2 and cm respectively. 
As a result the Feinstrukturkonstante a = hg 2 is dimensionless while £ 2 = hn is the Planck area. 
The kinematical phase space of the theory is subject to the (non - Abelean) Gauss constraints 

Gj = d a E a J + fj K L A^El 

Gj = d a E° + e jk l A k a E? (2.2) 

where in terms of a basis r j, Tj for Lie(G) and su(2) respectively the structure constants are de- 
fined by [zj,Zk] = fjK L t~l and [Tj,Tfc] = ej^ l T\ respectively. We normalize the anti Hermitean 
generators such that Tr (tjT k ) = —5jk/2 and Tr (t^t^) = —5jk/2. The phase space is subject to 
further spatial diffeomorphism and Hamiltonian constraints but these we reserve for a separate 
paper [To] . 

We now review the kinematical Hilbert space of the theory which is actually selected by 
requiring spatial diffeomorphism invariance It is easiest described in terms of so - called 
spin (or charge, flavour, colour..) network functions. We do this for a general compact group. 
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Definition 2.1. 

i) 

A spin network s for a gauge theory with compact gauge group G over a manifold a consists 
of a quadruple s = ( / y(s),Tr(s),m(s),n(s)) consisting of an oriented (piecewise analytic) graph 
7(5) embedded into a as well as a labelling of each of the edges e of 7(5) with a nontrivial 
irreducible representation ir e (s) of G and corresponding matrix element labels m e (s),n e (s) = 
1, .., dim(7r e (s)). Here for each equivalence class of irreducible representations we have chosen 
once and for all an arbitrary representative, 
ii) 

A spin network function is simply the following complex valued function on the space A of smooth 
connections over a 

T S (A):= [] [v^hnMi)){[7r e ( S )](^( e ))} me(sK(s) ] (2.3) 

where 

A(e) :=Ve*Veij A J tj) (2.4) 

denotes the holonomy of A along e and E{^{s)) is the set of edges ofj(s). 
Hi) 

The kinematical Hilbert space Hkiu is the closure of the finite linear span of spin network 
functions which define a basis. Hence the kinematical inner product is given by 

< T s ,T s i >Rin= $s,s> (2-5) 

The holonomy operators A{e) act by multiplication while the conjugate electric flux operators 
correspondinng to 

Ej(S) = [ (*Ej) (2.6) 
Js 

act by differentiation whose details JZ$ we will not need in what follows. Here *Ej denotes the 
pseudo two - form dual to the vector density Ej of weight one. 

One can show that TLKin = L>2(A, dfio) is a space of square integrable functions over a 
distributional extension A of A with respect to a probability Borel measure j!21 [H] but this 
will not be needed in what follows. For the Yang - Mills and gravitational sector respectively 
we have one kinematical Hilbert space each corresponding to the gauge group under question 
and SU(2) respectively and we will denote them as Ti^f^ and H-xfri respectively The total 
kinematical Hilbert space for the Einstein - Yang - Mills theory is simply the tensor product 
7~t-Kin = ^Kin ® ^xfrf ^ convenient basis is given by the states T c (S> T s where c, s respectively 
are spin networks for G, called colour networks in what follows, and for SU(2) respectively. 
Please refer to jS] and the sixth and eigth reference of [5] for more details. 

The functions (|2.2fl are scalar densities of weight one with respect to spatial diffeomorphisms 
Diff(cr) and they transform in the adjoint representation of G and SU(2) respectively. In [Hj 
we showed that their smeared form admits a well-defined quantization as essentially self- adjoint 
operators on TLxin- More specifically, let A, A be a Lie(G)— or su(2)— valued functions on a 
respectively (not necessarily smooth!) and let 

G(A) := / d 3 xA J Gj (2.7) 



4 



and similarly for G(A). Then 



G{A)T C ®T S = {ia Y, A J («)[ Y Stj- Yl ^ T c}®Ts 

v&V(-y(c) eeB(7(c));6(e)=-u e£E(~/(c));f(e)=v 

= ■ Y A J (^j T C }®T S 

v£V(l(c) 

G(A)T C ®T S = T c ®{i% Y, A >H E R J- J2 L ^ T ^ 

ueV(7(s)) eG£(7(s));6(e)=j; e£E(~f(s));f(e)=v 

=: T c ®{U 2 p Y A »*J T.} (2.8) 
uev(7(c)) 

Here V(7) denotes the set of vertices of a graph 7, 6(e), /(e) respectively denote beginning and 
final point of an edge e and the operators R%L e - respectively act as follows (the description 
is similar for the Yang - Mills counterparts): We may view spin network states as so - called 
cylindrical functions 

) (2-9) 
where / 7 ( s ) is a complex valued function on SU(2)' E ^ S ^. Then 

(R*T S )(A) := (^ = o/ 7(s) ({e^^A(e')}, eE(7(s)) ) 

mT s ){A) := (^)t=of l{ s){{A{e')e tT ^'}e^E { i{s))) ( 2 - 10 ) 

This ends our review of the kinematical LQG decription of Einstein - Yang - Mills theory. In the 
next two subsections we will construct and solve the associated Master Constraint Operator 
corresponding to both Gauss constraints. 



2.1 Einstein — Yang - Mills Gauss Constraint 

In j^l it was shown that background independent theories offer the possibility to define ultraviolet 
finite operators for classical integrals of scalar densities of weight one. The intuitive reason for 
this is as follows: If F = J d 3 xf(x) is the classical integral to be quantized then the density 
n scalar f(x) becomes an operator valued distribution of the form f(x) = ^2 a (S(x,x a )) n f a 
where the sum is over some index set depending on the theory and /. Here f a is an actual 
operator (not a distribution). Notice that the cassical density weight of / is correctly carried 
by the n — th power of the 5— distribution but unless n = 1 this is ill-defined and requires a 
point splitting regularization ("operator product expansion") with subsequent renormalization. 
This is actually the reason why in [2] we needed a trace class operator in order to split the 
points of the square (dE) 2 which is a scalar density of weight two. Splitting points requires a 
background metric which is not allowed in background independent theories, however, on the 
other hand in background independent theories the density weight comes out to be always unity 
because of spatial diffeomorphism invariance! Namely only integrals of scalar densities of 
weight one are spatially diffeomorphism invariant. Hence in background independent theories 
renormalization is not only not allowed but also not needed. The interested reader is referred 
to the second reference in j^j and to jS] for more details. 

In keeping with this spirit we want to define a Master Constraint which is the integral of 
a density of weight one, quadratic in the Gauss constraint and independent of a background 
metric. The simplest and most natural choice is 

M := f d 3 x ~ J ^f (2.11) 
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Here the spatial metric q a b is defined via its inverse det(q)q ab = E^E b 5^ k . Notice that the 
overall density weight of (|2.11|) equals unity indeed. In fact, (|2.11|) is both G— invariant and 
Diff(<r)— invariant. The choice (|2.11|) is natural because \J det(g) is the simplest scalar density 
of weight one which at least classically is nowhere vanishing. We could not have proceeded 
similarly with pure Maxwell theory because there the simplest scalar density is d ■ E itself 
which would have resulted in M = J d 3 x\d ■ E\. Not only would it be impossible to quantize 
this expression on Fock space, it is also not a classically differentiable function on the classical 
phase space and hence inacceptable. In non - Abelean Yang - Mills theory it is also possible 
to construct the "metric" Q a b whose inverse is defined by det(Q)Q ab = E a ; E b K 5 JK , however, 



the corresponding scalar density y / det(Q) which does not vanish identically if dim(G) > 3 is 
classically not constrained to be non - vanishing everywhere. Hence the fact that we have 
coupled gravity cannot be avoided. 

To quantize ()2.11[) we will proceed similarly as in 0|. Hence we will be brief, referring the 
interested reader to the literature. The idea is to make use of the identity 

[det(e)] 2 



1 



det(g) 



(2.12) 



where e a is the co - triad defined up to SU(2) transformations by q a b = 8j^e? a e^. Let V be a 
partition of a into mutually disjoint regions R. Then the integral 1)2.11(1 is the limit as V — > a 
of the corresponding Riemann sum 



M = 


iim y - Am 

V — >IJ ^ 

Rev 


GAR) = 


[ d 3 xGj(x) 
Jr 


e(R) = 


/ d 3 xdet(e)(x) 
Jr 


V{R) = 


/ d 3 xy / det{q)(x) 
Jr 



V(RY 



We now observe the further identity 

1 



hence 
e(R) 



e(R) = (-r / e jk i{A>(x),V(R)}A{A*(x),V(R)}A{A'(x),V(R)} 
K Jr 



(2.13) 



(2.14) 



V(mW " 0' J R e ^ Aj ^> V ( R ) 1/2 } A i Ak ^ V(R) 1/2 } A {A l (x), V(R) 1 / 2 } =: (^f~e{R) 



Thus (|2.13j) becomes 



M 



y 6 E S JK GAR)G K (R)e(Rf 



(2.15) 
(2.16) 



Rev 



The reason for writing the Master Constraint in this form is that all quantities involved in 
()2.16j) admit a well-defined quantization at finite partition V: First of all Gj(R) = G(A) with 
A K = XR^Ij where \R ls the characteristic function of R admits a well-defined quantization 
according to (|2.8[) . Next the volume V(R) of R is well-defined as a positive essentially self- 
adjoint operator on H-Kin, its explicit action on spin network functions being given by 



V{R)T S 



ve 



V( 7 («))nB \ 



1 

48 



E 

ei,e 2 ,e 3 GE(7(s));fe(ei)= 



ei kl e(e u e 2 ,e 3 )RfRl 2 R?\ 



■b(e2)=b(e 3 )=v 



(2.17) 
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Here we have assumed that all edges are outgoing from a vertex (split edges into two halves 
to do that) and e(ei, 62, e^) is the sign of the determinant of the matrix defined by the column 
vectors e\ (0), 62(0)63(0) in this sequence where ej(0) = b(ej), I = 1,2,3. 

It remains to quantize e(R) itself and to take the limit V — > a. To that effect, notice that the 
limit in (|2.13|) is independent of the choice of the sequence of partitions V. We may therefore 
without loss of generality assume that the partition is actually a triangulation consisting of 
tetrahedra R. In the limit V — > a each R can be described by a base point v(R) and a right 
oriented triple of edges ej(R) C dR, I = 1,2,3, incident at and outgoing from v(R). Then it is 
easy to see that in the limit V — > a the expression e(R) can be replaced by 

-8e'(R) = -86^6 J ^ i Tr(r,(^(e J (i?)))- 1 {^(ej(i?)),l/( J R) 1 / 2 }) x (2.18) 
xTr(r fc (A(e i ,( J R)))- 1 { J 4(eA'( J R))^(^) 1/2 }) Tr(r / (A(e L ( J R)))- 1 {^(eL( J R)), V(R) 1 / 2 }) 

Expression (|2.18l) can now be quantized by replacing Poisson brackets by commutators divided 
by ih. 

Anticipating that the strong limit V —* a exists as a symmetric operator we define 

M = lim {-f Y * JK ^)^) 2 GAR) (2-19) 

Notice that no ordering ambiguities arise because the gravitational and Yang - Mills degrees of 
freedom commute with each other. It will be sufficient to define Q2.18JI on the basis T C ®T S . 

In performing the limit a regularization ambiguity arises: The quantum limit when applied to 
T c <g> T s depends on the sequence V — ► a while the classical limit was independent of that choice. 
We will therefore proceed as in [H] and choose a sequence for each T c (g) T s individually which is 
justified by the fact that classically the choice of the sequence does not make any difference. The 
result is the following: We notice first of all from (|2.8|) that there is no contribution in (|2.18|) 
from those tetrahedra R which do not contain a ver texof 7(c). Next, [A(e), V(R) l l 2 } evidently 
vanishes for those R which do not contain a vertex of j(s) no matter how we choose e C dR. 
Hence, for sufficiently fine V we can focus attention on those R which contain a common vertex 
of j(s) and 7(c). The idea is now to average over all triangulations. Let n(v) be the gravitational 
valence of such a vertex then the detailed averaging performed in [S] results in the following final 
expression after taking the limit V — > a 

MT C ®T S = a 2 {jf b JK Xj{v)^y) 2 X K (v) T C ®T S 

p vev{y(c))nv(rf{a)) 

^0 = , w , n 1 w -e^ kl e JKL Y x 

KJ n(v)(n(v) -l)(n(v) -2) ^ 

ei,e2,e3G-E(7(s)); b(e\)=b(e2)=b{ez)=v 

xTr(r J (^( S (e / )))- 1 [^( S (e 7 )),y( W ) 1 / 2 ]) Tr(r fc (^( S (e^)))- 1 [A( S (eK)), V{vf/ 2 ]) x 
xTr(r / (A( S (e L )))- 1 [A( S (e L )),y(t;) 1 /2]) (2.20) 

We have displayed the action for vertices of gravitational valence at least three. For gravitation- 
ally bi-valent vertices the action looks similar and is non- vanishing These details will not be 
important for what follows. 

Here s(e) denotes an infinitesimal beginning segment of an edge e and V{v) denotes the 
volume operator for an infinitesimal region R(v) containing v. By the methods of |5] it is easy 
to see that (|2.2()[) is a positive, essentially self-adjoint, gauge invariant and spatially diffeomor- 
phism invariant operator on Ti-Kin- Moreover, the choice of the segments and regions s(e),R(v) 
respectively is completely irrelevant as long as they contain 6(e), v respectively. 
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To solve the Master Constraint (|2.2U|) is now surprisingly simple because the spectrum of 
the Master Constraint Operator (|2,2U|) is pure point. Moreover, we can determine it sufficiently 
explicitly in order to actually derive the physical Hilbert space. To see this, notice that "Hj^ 
can be decomposed as 

T^Kin = ©7^tn; 7 ,{7r} ( 2 - 21 ) 

where is the finite linear span of charge network states c with 7(c) = 7, tt(c) = 7r and 

the overline denotes closure. Now each Hj^.^ n can be further decomposed as 

^ifin; 7 ,7r = ©n^^;-,^^ (2.22) 

where II = {^lv} v ^v(j) is a collection of equivalence classes of irreducible representations of G, 
one for each vertex of v, with the following meaning: From IJ2.HJI one easily verifies that the 
colour network state T c , 7(c) = 7, 7r(c) = tt transforms under local G— gauge transformations 
with support at v £ V(j) in the representation 

[®ee£( 7 ); b(e)=v^e] ® [®ee£( 7 ); f{e)=vK] (2.23) 

where tt c denotes the representation contragredient to n. Since G is compact, every represen- 
tation is completely reducible and (|2.23f) can be decomposed into mutually orthogonal but not 
necessarily inequivalent irreducible representations. Let the orthogonal projector on the repre- 
sentation space consisting of mutually orthogonal representations equivalent to the equivalence 
class 11^ be denoted by 

l [®b(e)=*v*e]®[®f{e)=v*$ (2.24) 

also called an intertwiner. Notice that (|2.24|) vanishes for all n„ except for finitely many. We 
have the completeness relation 

^2 l f®b(e) = vKeM®f {e) = v lT<i] = ld [®i,(e)=^ e ]®[<g> /(e)=1) 7rg] (2-25) 

n 

Using (|2.25j) for every v £ we can decompose each colour network state and arrive at 
Finally we define 

C ;7 ,n = ©^™ ;7 ,n,. ( 2 - 26 ) 

and thus have the identity 

-nlfn = e 7l nW^. 7i n (2-27) 

The point of the decomposition Q2.27JI is that the operator A 7i „ := 5 JK XjX^ which appears in 
(|2.2(Jj) is diagonal on Hj^.^ n with eigenvalue — An„ < given by the eigenvalue of the Laplace 
operator on G in any representation equivalent to IT. This also demonstrates that the sum in 
(|2.27f) is indeed orthogonal. 

We can now proceed similarly with T~Lp^ n and decompose it as 



H%? n = ® 1 ,xH%f n . 7jA (2.28) 

where —A = {— A„}„ g y( 7 ) < are the eigenvalues of (e'{v)) 2 . That this operator is diagonizable 
in this fashion follows from a similar property for the volume operator itself. The direct integral 
decomposition of Tixin with respect to the Master Constraint Operator is therefore simply 



n Km = e (7 ,n), (7 ',A)^ ;7 ,n ® n%f n . j/>x (2.29) 

Since the spectrum of M is pure point we just need to identify the zero eigenspace in (|2.29|) as 
the physical Hilbert space which, as a subspace of H-Kim carries the kinematical inner product 
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as the physical inner product. To that end we notice that the eigenvalue on "H^^ n <8> 'H'xfn 'y \ 
is given by 

« 2 (^) 6 £ An„ (2.30) 
p tiev(7)nv(7 / ) 

If the > are not vanishing, the only way to make (|2.3()fl vanish for all 7 is to require An„ = 
for all v E V(7). Hence 11^ = Triv must be the equivalence class of the trivial representation in 
that case. 



If X v = then 11^ is arbitrary, which is different from what the ordinary Gauss constraint 
would select. We actually do not know explicitly the space of states with A„ = but it contains 
at least states with at least one vertex which is not at least trivalent with respect to the grav- 
itational spin network structure. These zero eigenvalues are related to zero volume eigenstates 
because e'(B) is related to the volume of the region B and we know by now that there are many 
such states even for higher valent vertices To see that this does not pose any problem, 

notice that a zero volume vertex physically corresponds to a region which actually does not 
exist. In other words, whether the Gauss constraint holds there or not is unaccessable by any 
observer 1 . Likewise, the results of |19j demonstrate that the number of zero volume eigenvalues 
is neglible in a semiclassical sense because they correspond to spin configurations whose number 
compared to the number of all spin configurations at an n— valent vertex with maximal spin 
j decreases roughly as j~^ n ~^ in the semiclassical limit of large quantum numbers j. Hence, 
semiclassical states are linear combinations of spin network states with almost always \ v > as 
it should be since classical General Relativity is about non - degenerate metrics. In what fol- 
lows we will simply drop that unphysical subsector of zero volume states for notational simplicity. 

Then, denoting 

1~Llnv := ©7^i^ ;7 ,Triv (2.31) 

the physical Hilbert space is given by the subspace 

Hphys = Klnv ® ^Kfn C ^Kin (2.32) 

consisting of the closed linear span of states T c ® T s where T c is a gauge invariant charge network 
state constructed by using the intertwiner for the trivial representation. This is precisely the 
same phyiscal Hilbert space as selected by the Gauss constraint (|2.8|) itself. 



2.2 Einstein Gauss Constraint 

We can proceed completely similarly with respect to the gravitational Gauss constraint for which 
the Master Constraint is given by 

M : = f d 3 x ^£= (2.33) 
Jo y/det(q) 

1 One could object that the zero volume region has a closed two - boundary for which we could construct the 
electric flux operator. However, notice that in non - Abelean gauge theories the electric flux operator for a closed 
two surface 1. does not measure the charge contained in the region bounded by the surface because it cannot be 
obtained by the (non - Abelean) Stokes theorem from the Gauss constraint and 2. it is not a Dirac observable 
even with respect to the Gauss constraint, that is, it is not gauge invariant and therefore does not correspond to 
a physical (measurable) quantity. The area operator is gauge inavariant but does not measure the charge neither, 
in fact, the charge is associated with vertices while the area is associated with edges. For Abelean gauge theories 
the objection is valid but here we are interested in a gravity coupled situation and then we have to consider also 
the spatial diffeomorphism and Hamiltonian constraints respectively for which the flux is not a Dirac observable 
again. 
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Its quantization proceeds entirely analogous to that of the Einstein - Yang - Mills Gauss con- 
straint, just that it acts trivially on the Yang - Mills sector so that we will drop 'hCjM, for the 
remainder of this section. Its action on spin network functions is given by 

MT S = ^(1) 6 J2 ^ k Xj{v)7{v) 2 X k {v) T s (2.34) 

Now one could imagine that an ordering issue arises: In the ordering (|2.34j) the operator is 
manifestly positive and essentially self-adjoint. However, in order to conclude similarly as in the 
previous section we need the ordering in which 8i k Xj(v)Xk(v) stands to the outmost right. In 
the Yang - Mills case there was no issue because the gravitational and Yang - Mills degrees of 
freedom commute. Fortunately also here there is no problem because Xj(v) is the generator of 

SU(2) gauge transformations at v and the operator e'(v) is manifestly gauge invariant. Hence 

-2 

the operators S 3 Xj(v)X}.(v) and e'(v) commute and can be diagonalized simultaneously. Hence 
we can perform exactly the same steps as in section |2~T1 and conclude that the physical Hilbert 
space is given by 

7~Lphys = 7~(-Inv C H-Kin (2.35) 

3 Infinite Number of Non-Abelean First Class Constraints Non 
— Polynomial in the Momenta with Structure Functions 

Euclidean 2+1 gravity can be formulated in complete analogy to Lorentzian 3+1 gravity (see e.g 
the fifth reference in [S]). In this form the constraints are as difficult to solve as for the 3+1 theory 
which we reserve for future work 15 . In particular, they are non-polynomial in the momenta 
and only close with structure functions rather than structure constants. However, there is a 
classically equivalent way (when the spatial metric is not degenerate) to write the constraints 
linear in the momenta. Namely Euclidean 2+1 gravity can be written as an SU (2)— gauge theory 
over a two dimensional Riemann surface a subject to a Gauss constraint and the curvature 
constraints 

& = \e ab Fi b (3.1) 

where F'K is the curvature of A 3 a and e ab is the density one valued skew tensor in two dimensions. 
It follows that (jHU is a scalar density of weight one. Accordingly we define the corresponding 
Master Constraint as 

ME : = l[ d ^E^M (3.2) 



2 J, v'dSh) 

where det(q)q ab = EjE^.5 jk defines the spatial metric. Of course the pair (A J a , Ej ) is canonically 
conjugate up to the gravitational coupling constant k. The significance of the label E will be 
explained momentarily. 

To quantize Q3.2JI we will work directly on the space of gauge invariant states Tli nv given by 
the closed linear span of gauge invariant spin network functions derived in the previous section. 
We proceed as in the fifth reference of jH] and define the "degeneracy vector" 

BP = \e ah ^ kl EtE\ (3.3) 

called this way because det(g) = 5 jk EiE k . Using this we can establish an interesting identity: 
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Notice that the combinations (we use the identity F 3 ab = e^F 3 ) 
C a := CjCabEj = Fl b Ej 

C l C 3 e ab EtE^ kl _ l F^EjEjeW 

y/det^J 2 yditM 2 ydetM 

look excacly like the spatial diffeomorphism and (Euclidean) Hamiltonian constraint of the 3+1 
theory. Now one immediately verifies that 

Me = vm (3 - 5> 

Hence we arrive at the following crucial observation: 

The Master Constraint (|3.2|) for the curvature constraints (|3.1|> coincides with the 
extended Master Constraint [lj for the spatial Diffeomorphism and Hamiltonian con- 
straints. 

Hence 2+1 gravity not only tests the extended Master Constraint idea but also can be consid- 
ered as an example with a non - Abelean constraint algebra with structure functions. 

Similar to the previous section we define for a two dimensional region R the following smeared 
quantities 



Fj (R) = [ d 2 xF, 

JR 

E° (R) = [ d 2 xE j 
JR 



V(R) = / d 2 x^det{q) (3.6) 



R 



and write the classical Master Constraint as the Riemann sum limit for a family of partitions 
V, that is 

_ ,. m 1 E ISffiW (3.7) 



Using the classical identity for 

1 

K 

we have 



E?{x) = ^e ab e jkl {A k a (x),V(R)}{A l b (x),V(R)} (3.8) 



vwi~ 2 = & hjki L {Ak ' v(R)1/4} A {Al ' v{R)1/4} =: & E i {R) 



(3.9) 



Hence the Master Constraint becomes 

lim 



M £ = ^imJ(i) 4 ^(F J ( J R)) 2 (^( J R)) 2 (3.10) 
Rev 



Again we specialize to a simplicial decomposition V . We single out a corner v(R) for each 
triangle R and denote by ej(R) the two edges of dR starting at v(R). Then 

M E = ^2(i)*J2&W)*@hW) 2 (3-11) 
^ CT Rev 

Fj(R) = Tr(TjA(dR)) 

jki sr x 



E'{R) = j kl Y, 



K,L=1,2 

xTr{T k A{e K {R)){A{e K {R)r\V{R) 1 ^}) TrM(e i (i?)){A(e I/ (i?))- 1 , V(R) 1 ^}) 
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Now (|3.11|) is written in terms of holonomies and the volume operator both of which admit 
well-defined quantizations on TLxin (see the fifth reference of [Hj for the 2+1 volume operator). 

There is a difference in formulating the Master Constraint Operator M for the Hamiltonian 
constraint and the extended Master Constraint Operator for the combined spatial diffeo- 
morphism and Hamiltonian constraint: As shown in Q] the former can only be defined, in its 
fundamental form, on the spatially diffeomorphism invariant Hilbert space H-Diff QUI °f solu- 
tions to the spatial diffeomorphism constraint while the second of course must be defined on the 
kinematical Hilbert space Tixin since it solves both types of constraints in one step. However, 
since M^; is nevertheless a spatially diffeomorphism invariant operator, it can be defined on 
Ti-Kin only if it does not change the graph of a spin network state on which it acts as was shown 
in detail QU]. We must keep this in mind when quantizing ()3.11j) . 

Our heuristic ansatz will be to define the following quadratic form in which we order the 
various operator factors judiciously (if it exists) 




(3.12) 



x < T s , Fj(R)E' m (R)T Sl > Km < T s ,,F k {R)E' n {R)T Sl > Kin where 

EyR)T Sl = ^ kl eKLyc 

K,L=l,2 

xT,{r k A{e K {R))[A{e K {R)y\V{R) l ' A \) TT{nA{e L {R))[A{e L (R))-\ V(R) 1/4 }) T ai 

and Fj(R) acts by multiplication. Notice that the sum over s\ is an uncountably infinite sum 
but we will see soon that for given s, s' only a finite number of terms contribute so that the 
interchange of the limit with the sum in the second step is indeed justified. One may wonder 
why we have introduced this insertion of unity at all. This is motivated by the theory in 3+1 
dimensions ^ EH where this insertion of unity is mandatory when working at the level of Ti. eh ff 
and where T s ,T s i are replaced by diffeomorphism invariant distributions. 

It remains to take the limit V — > cr. First of all, as the regions R shrink to points, every term 
in the sum over the R at given si vanishes anyway if 7(5) 7^ j(s') due to the orthogonality of 
the spin network functions so that both graphs must be subgraphs of j(si) D dR, so they must 
equal each other for sufficiently small R. Thus: 

The M^; regularized this way in terms of a quadratic form is automatically not 
graph changing. 

Next, no matter how the limit V — > cr is performed, given s, for sufficiently fine V the last 
line in (|3. 12f) vanishes unless R intersects 7(5). Now in contrast to the situation in 3+1 di- 
mensions, if R intersects 7(5) but does not contain a vertex of 7(5), the last line in (|3,12|) is 
not automatically vanishing which leads to a divergence when the limit V —* cr is not carefully 
performed. Following the proposal in the fifth reference of [2j we adapt the limit V — » a for the 
matrix element ()3.12|) to the states T S1 in the following sense: 
1. 

The triangles R are chosen to saturate each vertex v of 7(si), that is, for each adjacent pair of 

edges e, e' there is precisely one triangle R with v(R) = v and ei(R) C e, e2(R) C e'. 

2. 
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Away from the vertices, the triangles R intersect 7(5) only in such a way that the vertices of 
the graph dR n 7(5*) which are not vertices of 7(s) are co-linear, that is, the tangents of the 
respective intersecting edges are linearly dependent there. 

This results in a different limit V — > a for each 7(si) which is justified by the fact that classi- 
cally the Riemann sum converges to the same limit no matter which sequence is chosen. The 
second requirement now makes sure that the only contributions come from the triangles with 
v(R) € V(7(s)), hence (|3.12j) simplifies to 

Qm e (T s , T 8 t) = 5 7(s)i7(s0 ]Thm2(l) 4 £ £ frx 

«1 a p d€V(7(s)) R&V; v(R)=v 

x < T S ,F^)E!JR)T S1 > Kin < T s ,,F k (R)E> n {R)T Sl > Kin (3.13) 

Now a miracle happens: The limit involved in (|3.13j) is already trivial because, as the trian- 
gles R(v) shrink, this shrinking can be absorbed by an (analytical) diffeomorphism of a which 
preserves 7(5) = j(s'). However, since the diffeomorphism group is represented unitarily on 
7~(-Kin jlQJ, the number on the right hand side of (|3.13|) actually no longer depends on the "size 
and shape" of the triangles R but only on their diffeomorphism invariant characteristics (this 
is easiest to see by performing the sum over s\ thus resulting in a unit operator). This reduces 
the number of regularization ambiguities tremendously from uncountably infinite to countably 
infinite (namely the C n class of the intersections at the two vertices of dR H 7(s) different from 
v{R)). 

Hence, we are left with making that choice. The most natural choice is the one made in 
[S] namely such that the the edge dR — 7(si) intersects 7(^1) transversally at both end points 
(corresponding to a C° intersection). However, it is easy to see that the limit in (|3.13j) is not 

satisfactory unless dR C j(s). The reason is that up to the commutator between F)~(R) , E' n (R) 
(which is of higher order in K) after summation over s\ the matrix element is proportional to a 
sum of matrix elements between spin network states over 7(5) of the operator 

(F^R))^F^R) = 3 - XliMdR)) 
where Xj is the character of the irreducible representation of SU(2) with weight j. Since 7(5) = 
j(s') and since E'-{R) only acts on the intertwiners of a spin-network state but not the spins, it 

follows that the piece xi(A(dR)) drops out of (|3.13|) for sufficiently small dR. Hence the matrix 
element would not contain any information about the curvature Fj. This happens because we 
are working at the kinematical level and not at the spatially diffeomorphism invariant level: If 
T s would be diffeomorphism invariant and depend only on the diffeomorphism class of s ^5] 
then there would be a non-vanishing contribution no matter how small dR because we could for 
instance choose s to be in the class of the spin network state xi(A(dR))T s i . 

It follows that at the level of TLxin we must proceed differently than for Ti-Diff in order to 
obtain a satisfactory operator. The necessary idea was outlined already in [Tl I16j: 

Definition 3.1. 

Let 7 be a graph, v £ V{l) 0, vertex and e,e' € E(j) two different edges starting at v (reverse 
orientation if necessary). A loop a 7j „ ieje ' within 7 starting at v along e and ending at v along 
(e') _1 is said to be minimal if there is no other loop within 7 with the same properties and less 
edges of '7 traversed, 
it) 

Given the data 7, v, e, e' we denote by L(j, v, e, e') the set of minimal loops compatible with those 
data. 
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Notice that the notion of a minimal loop does not refer to a background structure such as a 
metric. The set L(j,v,e,e') is never empty for the closed graphs that we are considering here 
and it is always finite. 

The idea is now that semiclassical states must necessarily have very complex graphs in order 
that, e.g. the volume operator for every macroscopic region has non-zero expection values. 
Moreover, they will involve very high spins because high spin means large volume so that the 
correspondence principle is satisfied. For such complex graphs the limiting loops dR are of the 
same order of magnitude as the minimal loops associated with it. Hence on semiclassical states 
it makes sense to replace the limit by those minimal loops. This is justified in the sense that 
we are quantizing an operator corresponding to a classical quantity which we can hope to be 
approximated well for semiclassical states only. Since then there is no consistency check for the 
details at the microscopic level, our procedure is justified in this semiclassical sense. Of course, 
the argument is not entirely satisfactory and can at best result in an effective description M E, e f / 
of the fundamental operator which, however, must then be defined on H-Diff rather than 
H-Kin in order to obtain a non-trivial result. Hence, at the level of Ttxin this is the best we can 
do. The more fundamental programme will be carried out in [T5: directly for the 3+1 case. 

Summarizing, the close to final matrix element is defined by 

QM E ,eff( T s,T s i) = 5 7 (s), 7 (s') x (3.14) 

*E 2 4) 4 E ^ mn E E 



11 £^ ^ ^ ^ I L( 'yi s) v Ci 69 ) I 

si P »eV( 7 (s)) ei,e 2 eE(7(s));u=eine 2 a£L(-y(s),v, ei ,e 2 ) 1 ww ' ' x ' 



x < T s ,F J (R(a))E' m {v,e l ,e 2 )T Sl > Kin < Tj, F k {R(a))E' n { V, ei, & 2 )T Sl > Kin 
E'^Xe')T Sl = e^ M £ e KL x 

K,L=\,2 

xTr(r^( S (e i ,))[A( S (e^))- 1 ,F(i;) 1 / 4 ]) Tr{ Tl A{s{e L ))[A{s{e L ))-\ V(u) 1 / 4 ]) T s 

Here we have averaged over the number of minimal loops with data 7(s), v, e, e' in order not to 
overcount as compared to (|3.13|) and R(a) is the unique, interior, contractable region enclosed 
by q, i.e. dR(a) = a. As before s(e) denotes an infinitesimal beginning segment of an edge 
e and V{v ) the volume operator for an infinitesimal region containing v. Notice that H3,14|) is 
a drastic modification of (|3.13|) unless the graph 'j(s) is very complex, hence we have put the 
subscript "effective" . 

Formula (|3.14|) is not yet quite what we want because we have to impose the constraint 
&y( s ),7(s') explicitly. We can avoid that by inserting a projection operator P a which projects 
onto the closed subspace of spin network states which have spin higher than 1/2 for all edges 
of a. This again modifies M^for spin network states involving low spins but does not change 
its semiclassical behaviour for the reason mentioned above. Hence the final expression for the 
effective Master Constraint Operator is 

Q MB , e// (T s ,T s ,) = J>(1) 4 E ^ m " x ( 3 - 15 ) 

si P ueV(7(s)) 



E E 



\L(j(s),v, ei , 62)1 

ei,e26-Z?(7(s));u=eine2 a£L(-y(s),v,ei,e2) 



x < T s ,F j (R(a))E' m (v,e 1 ,e 2 )P a T Sl > Km < T^F k (R(a))E' n (v,e 1 ,e 2 )P a T 8l > Ki 
Due to the projections the state 

F j {R(a))E' m (v,e 1 ,e 2 )P a T Sl 
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is a linear combination of spin network states over the graph 7(«i) U a = 7(^1 ), hence we must 
have 7(5) = 7(si) = j(s') in order that the matrix element for s\,a does not vanish. Thus the 
condition 7(5) = 7(s') is implicit. 

Expression (|3.15j) defines a positive quadratic form. To see the positivity we take a generic 
linear combination of spin network functions 

M N p 
p=l r=l 

where z r are complex numbers and we have adopted a labelling such that 7(sp = 7 P for all 
r = 1, ..,iVp and the graphs 7 P are mutually different. Then we simply compute, using the fact 
that the matrix elements vanish between spin network states over different graphs 

N N p 

Qm b , £// (/,/) = E E "p z p x ( 3 - 16 ) 

p=l r,s=l 

*E 2 <|> 4 E E E y 



x < T s r,F i (i?(a))£ , ^(i;,ei,e2)-P a T Sl >^ in < TU, F k (R(a))E' n (v , ei, e 2 )P a T Sl > Ki 



N 



EE 2 9 4 EE E E 

\L(j p , v, e\, e 2 )\ 



■:iE^ < T s r,F j (R(a))E^ n (v,e 1 ,e 2 )P a T Sl > Kin \ 2 



r=l 

which is manifestly non - negative. 

Not every positive quadratic form defines an operator but if it does (technically, if it is 
closable), then the corresponding positive, self-adjoint operator is unique (i.e. there is no choice 
in its self-adjoint extension). 

Theorem 3.1. 

The quadratic from A ff. 15)1 defines a positive self - adjoint operator. 

Proof of theorem 13.11 
To see this we notice that the would be operator is given by 

M E ,effT s = J2QM E ,eff( T s',T s )T s > (3.17) 

a' 

This defines a (densely defined) operator if and only if the right hand side of (|3.17ft is normal- 
izable, that is, 

\\ME,effT s \\ 2 Kin = Y,\QM E , eff (T s >,T s )\ 2 (3.18) 
s> 

converges. Notice that the sum on the right hand side is over an uncountably infinite set, so 
()3.17|) looks dangerous. However, we know already that ()3.18p reduces to the countable sum 

II ^E,effT 8 \\jcin = E \QM Eteff {T s r,T s )\ 2 (3.19) 
7(s')= 7 (s) 

Now the matrix elements (|3,15[) are finite for given s,s' since we must have 7(^1) = 7(5) and 
the states Fj(R(a))E' m (v,ei,e 2 )P a T Sl are a finite linear combination of spin network states of 
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which at most one can coincide with either T s or T s i. Since there are only Y^veVMs)) n ( v ) such 
terms involved in (|3.15|) where n(v) denotes the valence of v, the assertion follows. 

Hence, to show convergence of (j3.19|) it is sufficient to show that for given s there are only a 
finite number of s' for which Qm e eff (T s i ,T S ) ^ 0. Now we have seen that Qm e eff (T s i ,T S ) ^ 
implies that in the sum over s\ only a finite number of terms contribute which must satisfy 
s\ € S(s) where S(s) is a finite set of spin network labels. For each s\ £ S(s) consider the finite 
number of spin network states with labels in S(si) contained in the span of the states 

Fj(R(a))E' m (v, e 1 ,e 2 )P a T Sl 

as v,e\,e2,a € L{ r y(s), v, ei, 62) vary. Thus s' must be in the finite set U 81 ^gf 8 \S(si). 

This shows that (|3.17|) defines a positive symmetric operator on the finite linear span of spin 
network states. However, positive symmetric operators have a preferred self - adjoint extension, 
the Friedrichs extension, which is given by the unique closure of the corresponding quadratic 
form Qm e eff with the same domain of definition. 
□ 

We now can proceed to solve M Et eff- Let &Kin be the dense subset of Hxin given by the 
finite linear span of spin network functions. We begin by verifying that for any / G &Kin the 
distribution 

n(f):= fS[F], S[F}:= f du (A ) 5 Aq (3.20) 

J M flat 

solves the effective extended Master Constraint Operator. Here M. fi a t is the moduli space of 
flat connections on a, 5a is the 8— distribution supported at Aq and v$ is the following measure 
on Mfiat- Any function on Mfi at is of the form f(A Q ) = f n (A (ai), ..,a n (A )) where a l7 ..,a n 
are any generators of the fundamental group vri(cr). Then 

Mf)= [ dis (A )f(A )= f dii H {hi)..dix H {K) f n {hi,..,K) (3.21) 

J M nat JSU(2)« 

The interested reader is referred to [S] for more details. 

Now we may check whether (|3,20[) is a generalized solution of ME,eff = 0. We have for any 



v(f)[M E>eff f] = [ du (Ao)f(A ) S Ao [M E>eff f) 

J M f l at 

= [ dv (A )J(Aoj[M Eteff f'](A ) (3.22) 
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This will be zero if and only if it is zero for any spin network state /' = T S , hence we compute 

V(f)[M E ,effTs} = Yl < f> T s> >Flat Q& E „(?>,T S ) 



p s',si DeV(7(s)) 

E E 



iLfjfs), v, e%, e 2 ) 

e\,e2EE(y(s));v=ei_r\e2 ae.L(-y(s),v,ei,e2) 



x < T s i ,F j (R(a))E! m (v,e 1 ,e2)P a T Sl > Kin < T s ,F k (R(a))E^(v,ei,e 2 )P a T Sl > K i 

2(1) 4 V V P k 5 mn V V — --i ^x 

B ^ ^ t-^i \Lh(s),v,ei,e 2 )\ 

F s i v£V(-y(s)) ei,e2€^(7(s));«=eine2 a6L(7(s),«,ei,e2) 



x < f,F j (R(a))E' m (v,e 1 ,e 2 )P a T Sl > F iat < T s ,F k (R{a))E' n {v,e 1 ,e 2 )P a T Sl > K in 
= (3.23) 

In the last step we have used the completeness relation with respect to s' and the fact that 
Tr(TjAo(a)) = for every Ao € M.fi a t- We could interchange the various summations because 
the non-vanishing terms reduce the sums to finite ones. 
It follows that 

M' E , eff v(f) = (3-24) 

where the prime denotes the dual of the operator on the space 3>|f » n of linear functionals on &Kin 
(without continuity requirement) defined in general by [0'l](f) = ZfO'/]. Hence we see that the 
point A = lies at least in the continuous part of ME,eff since the rj(f) are not normalizable 
with respect to < ., . >Kin, provided that we can write ??(/)[/'] =< /(0), /'(0) > w =e where /, /' 
are representatives of /, /' corresponding to a direct integral decomposition of TL subordinate to 
M^, One should now complete the analysis and compute the full spectrum of Mg^ff to see 
whether that is indeed the case. Assuming that to be true, for the purposes of this paper it is 
sufficient to note that whatever the complete Tiphys might be, it contains the closed linear span 
of the r](f) as a subspace with the induced physical inner product 

<v(f),v(f')>Phy S :=Mf7) (3.25) 

which is of course well known in the literature |20| . We leave a more complete analysis for 
future work and just remark that with the techniques of JZ] it is not difficult to show that 
the operator (|3.17|) . although we have performed rather drastic manipulations, indeed has the 
correct classical limit. 



4 Conclusions and Outlook 

What we have learnt in this paper is that the Master Constraint Programme is also able to deal 
with the case of interacting quantum field theories. By this we mean that we can solve, e.g. the 
Gauss contraint of Non - Abelean Yang Mills theory when coupled to gravity. Classical Non - 
Abelean Yang - Mills theory on a background spacetime is a self - interacting field theory and 
in four dimensions nobody was able to show that the corresponding interacting quantum field 
theory exists (in the continuum; on the lattice there are no problems). The Gauss constraint of 
the theory is a quadratic polynomial in the fields and its square id a fourth order polynomial. 
Therefore, from a QFT on curved background point of view the Master Constraint should be as 
UV singular as the square of the curvature that defines the Yang - Mills action. Things become 
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even worse because on top of that we have multiplied the square of the Gauss constraint by 
a factor which depends non - polynomially on the the degrees of freedom of the gravitational 
field. 

Yet, we were able to quantize the integral of the resulting expression on the Hilbert space 
which is used in Loop Quantum Gravity (LQG) without encountering UV problems. The tech- 
nical reason for why this happened is that the Master Constraint is spatially diffeomorphism 
invariant integral of a scalar density. As was shown in the sixth reference of for such quan- 
tities, in a very precise sense, the UV regulator gets swallowed up by spatial diffeomorphism 
group: The theory does not depend on a background metric, thus all "distances" are gauge 
equivalent. 

Let us compare the situation with our previous paper [5]: There we were looking at free field 
theories on a Minkowski background and we could use the associated background dependent 
Fock representations. The square of the Maxwell Gauss constraint is too singular and cannot be 
employed to define the Master Constraint, thus we had to use the flexibility of the Master Con- 
straint Programme to regularize the Master Constraint by a background dependent kernel. In 
the present situation the singularity of the square is expected, according to perturbation theory 
arguments, to be even worse. However, fundamentally the Fock representation is not a valid rep- 
resentation for interacting quantum field theories. A possible representation is the one that uses 
LQG techniques, however, that representation is only valid when we couple the gravitational 
field as otherwise e.g. the Yang - Mills Hamiltonian is ill - defined in that representation. 

It is precisely this observation which has lead us to quantize the Master Constraint in this 
kind of LQG representation. It turns out that the non - polynomial, gravitational field dependent 
function mentioned above that enters the Master Constraint plays exactly the same role as for 
the Yang - Mills Hamiltonian (constraint): It serves as a background independent UV regulator. 
Moreover, while in the free field case the regulator could be chosen to be background dependent, 
here we are not allowed to do that, however, background independent theories, through gravity, 
have the tendency to regulate themselves. Thus the factor of 1/ \J det(q) which typically enters 
the background independent Master Constraints becomes a quantum operator (rather than a C 
number valued expression) but otherwise plays the same role as the trace class operator K in 
[3]: It removes the UV singularities of the square of the Gauss constraints. 

We stress again that this does not mean at all that we have proved the existence of, say, 
QCD. This is because, while we are able to avoid certain singularities, now the burden is on us to 
show that the theory can also successfully deal with the additional symmetries that have entered 
the stage by coupling matter to the gravitational field. This is the spacetime diffeomorphism 
symmetry which finds its way into the canonical framework in the form of the spatial diffeomor- 
phism and the Hamiltonian constraint. It is precisely for this reason that the Master Constraint 
Programme was created. One now has to apply it to all symmetries of General Relativity, solve 
the full Master Constraint and establish that we have captured a quantum theory of General 
Relativity rather than a mathematically consistent but physically uninteresting quantum theory 
of geometry and matter. This is what has to be done in the close future and finally the mathe- 
matical techniques are available in order to make progress. 
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